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Abstract 

We studied the transport properties of electrons in graphene as they are scattered by a double 
barrier potential in the presence of an inhomogeneous magnetic field. We computed the transmission 
coefficient and Goos-Hanchen like shifts for our system and noticed that transmission is not allowed 
for certain range of energies. In particular, we found that, in contrast to the electrostatic barriers, 
the magnetic barriers are able to confine Dirac fermions. We also established some correlation 
between the electronic transmission properties of Dirac fermions with the Goos-Hanchen like shifts, 
as reflected in the numerical data. 
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1 Introduction 


Graphene, a planar arrangement of carbon atoms on a honeycomb lattice, is a unique realization of 
a two dimensional electronic system. Due to its excellent carrier transport properties, graphene has 
a great potential for nano-electronic applications. Among the peculiar electronic properties of this 
2D-material is its unusual quantum Hall effect [1]. Graphene is also a transparent conductor [2] whose 
carriers are massless and chiral relativistic fermions governed by a Dirac-like equation leading to many 
fascinating physical properties of graphene, such as Klein tunneling [3,4]. However, as appealing as 
the Klein tunneling may sound from the fundamental research point of view, its presence in graphene 
is unwanted when it comes to applications of graphene because space confinements of the carriers 
is of great importance in nanoelectronic applications. In addition, the ability to control electronic 
properties of a material by an externally applied voltage is at the heart of modern electronics [5,6]. 

The inability to confine electrons using an electrostatic potential barrier severely limited the appli¬ 
cability of graphene based devices. However, it came as a big relief when it was pointed out that well 
localized magnetic field dubbed as magnetic barrier can confine massless Dirac fermions in graphene [7]. 
Later on, snake states, trajectories of charge carriers curving back and forth along interfaces, were 
proven to play an important role and were studied experimentally [8,9], mainly motivated by the quest 
for electrical rectification. The inhomogeneous magnetic field case in graphene was analyzed in [7]. 
Theoretically, electron waveguides, in graphene subject to a suitable inhomogeneous magnetic field, 
were considered in [10]. One of the interesting features of such inhomogeneous magnetic field profile is 
that it can bind electrons, contrary to the usual potential step. Such a step magnetic field will indeed 
result in electron states that are bound to the Hj-field step and are able to move only in one direction, 
along the step. 

During the past few years there was substantial progress in studying electron transport properties 
in graphene, among these developments we cite the quantum version of the Goos-Hanchen effect 
originating from the reflection of particles from interfaces. Many works on various graphene-based 
nanostructures, including single barrier [11], double barrier [12,13] and superlattices [15], showed that 
the Goos-Hanchen like (GHL) shifts can be enhanced by the transmission resonances and controlled 
by varying the electrostatic potential and induced gap [11]. Similar to the situation in semiconductors, 
the GHL shifts in graphene can also be modulated by electric and magnetic barriers [16], and atomic 
optics [17]. It has been reported that the GHL shifts have a major effect on the group velocity of 
quasiparticles along interfaces of graphene p-n junctions [18,19]. 

Very recently, the GHL shifts for Dirac fermions in graphene scattered by double barrier structures 
have been studied in [13]. Moreover, in [14] we have explored the zero, positive and negative quantum 
GHL shifts of the transmitted Dirac carriers in graphene through a potential barrier with vertical 
magnetic field. Numerical results show that only one energy position at the zero GHL shift exists 
and is highly dependent on the y-directional wave vector, the energy gap, the magnetic field and the 
potential. The positive and negative GHL shifts happen when the incident energy is more and less 
than the energy position at the zero GHL shift, respectively. In addition, we found that there are 
two values of potential at the zero GHL shifts, where a potential window can always keep the positive 
GHL shifts. These results may be useful in designing a graphene-based valley or spin splitter as well 
as manipulating the electrons and holes in graphene nanostructure. 
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Motivated by different developments on the subject and in particular as a follow up on our recent 
works [13,14], we investigate the GHL shifts in a gaped graphene system in the presence of an inho¬ 
mogeneous magnetic field and a double barrier potential. We separate our system into three regions 
and determine the solutions of the energy spectrum in each region. Matching the wave functions at 
both interfaces, we then calculate the transmission coefficient as well as the GHL shifts. To allow a 
better understanding of our results, we study the transmission coefficient as well as the GHL shifts 
while varying different physical parameters that characterize our system. 

This paper is organized as follows. In section 2, we formulate our system Hamiltonian describing 
particles scattered in graphene by a double barrier potential in the presence of an inhomogeneous 
magnetic field. We then obtain the solutions of the energy spectrum corresponding to each region 
in terms of different physical parameters and analyze the energy conservation law. In section 3, the 
scattering problem for Dirac fermions will be solved using continuity at the boundary, which will 
enable us to calculate the transmission coefficient and corresponding phase. The condition for full 
reflection are then obtained for certain incidence angles In section 4, we study the GHL shifts and 
transmission coefficient as well as discuss our main results. We present our main conclusions in the 
final section. 


2 Theoretical model 


We consider a system of massless Dirac fermions moving through a strip of graphene and subject to 
a potential, which has the form shown in the Figure I. The system contains five regions denoted by 
the index j = 1,2-•• ,5. The left region (j = I) describes the incident electron beam with energy 
E = vpe and incident angle 4>i where vp is the Fermi velocity. The far right region (j = 5) describes the 
transmitted electron beam with a lateral shift St and angle (j)^ but in the presence of an inhomogeneous 
magnetic field. We introduce in the intermediate regions j = 2, A and middle region j = 3 two different 
magnetic fields B 2 and H 3 , respectively, such as 


Bj{x) 


B 2 , di <\ X \< d2 
< B3, I X |< di 

0 , otherwise. 

\ ’ 


( 1 ) 


In the present study, we consider the system in an inhomogeneous magnetic field given by the config¬ 
uration (I) in addition to the presence of an energy gap tj in the regions 2, 3 and 4 defined by 


t2, dl <1 X |< ^2 

< t^, I X |< di 

0 , otherwise. 

\ ■ 


( 2 ) 


In order to study the scattering of Dirac fermions in graphene by the above double barrier structure 
we first choose the following potential configuration 


di <1 X |< d2 
I X |< di 
otherwise 


Vj{x) = 


V 2 , 

V3, 

0 , 


2 


( 3 ) 






Figure 1: Schematic diagram for Dirac fermions in an inhomogeneous magnetic field and passing through 
a graphene double barrier, with height V2 in the region di <| x |< ^2 and height ^3 in the region | x \< di. 
(a) the dashed lines show smooth electric potentials having error function distributions, (b) describes the 
incident, reflected, and transmitted electron beams with a lateral shift St- 


where j labels the five regions indicated schematically in Figure 1 that shows the space configuration 
of the potential profile. The Hamiltonian for one-pseudospin component in the j-th region can be 
written as 

Hj = Vpa • TT -|- Vj{x)l2 -h (^2 “ (4) 

where 0 is the Heaviside step function, tt = p-\- eAj /c is the two-component kinetic momentum with 
the canonical momentum p = —ih{dx,dy)'^, cr = {ax, cry) and az are the nsual Panli matrices, I 2 is 
the 2x2 unit matrix. Choosing the Landan gauge we select the vector potential A = (0, Ay, 0)^ 
that creates the inhomogeneons magnetic field defined by (1), imposing the continuity of this vector 
potential at the boundaries of each region requires that 

X < —d 2 
— d 2 < X < —di 

\ x\<di (5) 

di < X < d 2 
X > d 2 


Ay{x) = Aj = - X < 


ih{di - d2) - T^di, 

^2 B 3 

^x + (^ - ^)di, 
‘S 2 B2 S 3 


^X - - ^)di, 

‘sa Sj ‘S 3 

I^{d2 - di) + ^di, 

sa ‘B 3 


where the local magnetic length is defined hy Ipj = -s/cjeBj in our selected system of units {h= 1). 
The eigenvalues and eigenspinors of Hj in regions 1 and 5 are generated by the Dirac Hamiltonian 




0 VF [pxj 

VF [Pxj +i{py + lAj)] 


i [py + f Aj)] 
0 


( 6 ) 


and the time independent Dirac equation for the spinor il>j{x,y) = {(p^,(pj associated with energy 
E = vpc is given by 




which can be written as two linear differential equations of the form 



)=d A 

J 

' V fj 


(7) 


Pxj -i[Py + -Aj 


Pxj +i[Py + -Aj 


if. =eip 


= ^^j 


3 


( 8 ) 

(9) 


























The corresponding energy eigenvalues read as 


where the symbol sj = sign(e) and 


e = + {Py + 

Pxj = ^je^- (py + • 


( 10 ) 


( 11 ) 


with incoming momentum pj = {pxj,Py) and position r = {x,y). The incoming wave function takes 
the form 

/ \ 

1 


“ ^/2 ( 


^ipyr 


■'Pxj 


and Zp^. is given by 


^Pxj — 


Pxj + i {Py + f. -^j ) 


= SjC 


( 12 ) 


(13) 


\/ {PxjY + {Py + l^i)^ 

/ p —^A \ 

where sq = sgn(e) and (pj = arctan f ° ^ j is the angle that the incident electrons make with the 
x-direction, pxi and py are the x and y-components of the electron wave vector, respectively. The 
eigenspinors read as 


V2 ( z, 


JiPxjX+PyV) 


(14) 


V2 ( -4 


J{-PxjX+Pyy) 


(15) 


To be much more accurate, we give the solutions of the energy spectrum for each region. Then in 
region 1 ( x < — ^2 ), we have 


e = 


pIi + 




Py + d 2 ) - -j^di 

^B2 ‘-Bs 


n 2 


1 / 1 


^1 = e^{Pxix+Pyy) ^ ^ 


1 / 1 


Ji-pixX+Pyy) 


V2\zi 

Pxi + i 


\/2 I -zf 


Zl = Si- 


Py + d 2 ) - ^di 

B2 B 3 


Iph + \py + J^idi - d 2 ) - J^di 

B 2 S 3 


and in region 5 (x > ^ 2 ); the solution is 


e = 


\ 


Px^ + 


Py 7^(*^2 ~ di) + j^di 


1 2 


T 1 / 1 

^5 = 

\/2 y Z 5 

Px5 + i 
Z5 = S 5 , 


^B2 


J{Px5X+Pyy) 


^B3 


Py + l^{d 2 - di) + 1 di 


l2 


Iph + \py + - di) + ^di 

B 2 


S 3 


(16) 

(17) 

(18) 


(19) 

( 20 ) 

( 21 ) 
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For the system under consideration, we can write the Hamiltonian corresponding to regions (2), 
(3) and (4) in matrix form as 


/ 


Hj = vp 


Vp ' Vp 


—I 




Ib^ 


J 2 c^-j) 




H i 

■^F "^F 


( 22 ) 


Note that the energy gap tj is equivalent to a mass term, this will lead to interesting consequences 
on the physical properties of such system. We determine the eigenvalues and eigenspinors of the 
corresponding Hamiltonian H by solving the time independent equation for the spinor ipj{x,y) = 
. Since the transverse momentum py is conserved, we can then write the wave function as 
'ipj{x,y) = e^yyytpj{x), with (pj{x) = and energy E = which lead to 


(23) 


At this stage, it is convenient to introduce the concepts of annihilation and creation operators in order 
to ease the diagonalization of our Hamiltonian. They can be defined by 

{dxj + ky+ = —k + ky+ 



)=d / 

^3 J 

' V 


V 2 


and obey the canonical commutation relations 


(24) 


T 


= 6 j^k- Rescaling our energies t- = vpPj and 


Vj = vpVj, then (23) can be written in terms of aj and a| as 


Vj + pj 

3 

^3 - H 




3 = 


= e 


^3 


(25) 


giving rise to the two relations between spinor components 

'Ib 


a /2 

{vj + = e(p+ 




+ (vj - = etp- . 


(26) 

(27) 


Injecting (27) in (26), we obtain a second order differential equation for 


/ \2 2 
{e-Vj) - y. 


+ ^ t + 


s, 


(28) 


which shows clearly that ip'J is an eigenstate of the number operator Nj = ajoj and therefore we 
identify (p'^ to be eigenstates of the harmonic oscillator \nj — 1), namely 


99 + ~| rij - 1 ) 


which is equivalent to stating 


/ \2 2 

[e-Vj) -//,■ 


^3 - 1 ) = I ^3 - 1 ) 


(29) 


(30) 


5 













and the energy spectrum can be defined by 


e-Vj — sjenj 


i = V 


where we have set = Sj (e — vj) and Sj = sign — vj) corresponding to positive and negative 
energy solutions. The second spinor component now reads as 






After normalization we arrive at the expression for the positive and negative energy eigenstates 


— - ^—7 -^ rtj — 1 ) 

^rijlBj ' -J ' 


\f 2 . lenjlBj-SjfljlBj 




CE^ / \ 

Introducing the parabolic cylinder functions Dn^{x) = 2“^ e~^Hnj ( j to express the solution in 
regions 2, 3 and 4 as 




/ ^rij ^Bj +Sj ^j^Bj 


iiSo \/2 






. (* 


TT-h kylSi 


FT 

3 


Ig. “*■ 


In summary the solutions of the energy spectrum in the barrier (—^2 ^ x < —di) (region 2) are 

€2 = U2 + S2^\l {y2lB2Y + 2^2 

^B2 


'Ip 2 {x,y) = a2i^t + ^ 2 V '2 


while in region 3 (|x| < di) read as 


€3 = ^3 +S3-— \nyslBs)'^+ ‘^n3 (37) 

'>p3{x, y) = + ftsV'g (38) 

and finally in region 4 (di < x < ^ 2 ) it can be expressed as 

€4 = V 2 + 547 ^ a/ {y2lB2Y + 2^4 (39) 

^B2 

'ipi{x, y) = + b4i>l (40) 

where the parameters aj and bj, with (j = 2,3,4), are normalization constants. 

Recall that, from the above analysis, we ended up with different energy spectra 62 , 63 and 64 , 
which are obtained in terms of system parameters and quantum numbers in each regions. On the 
other hand, energy conservation requires that 


e = £2 = £3 = €4 
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and by replacing the energies by their expressions, it is easy to observe that the allowed energy values 
should satisfy the relation 


77,2 = 774 


/2 

*52 


773 - 772 + + 


2773 




(42) 


Having obtained all solutions of the energy spectrum, we will see how they can be used to investigate 
different physical properties of our system. Specifically, we evaluate the transmission and reflection 
amplitudes in terms of different physical system parameters. 


3 Transmission and phase shift 


Before determining explicitly the transmission coefficient and its associated phase shift, we notice that 
total internal reflection will take place only when 0 < ijii < ^, since the wave incident from the right- 
hand and left-hand side of the normal surface will behave differently [20]. It is clear that the shift in 
Py is due to our choice of gauge for the vector potential. We find it more convenient to parameterize 
the momenta by 


Pxi = ecos(/>i, 
Pxb — €C 0 S(/) 5 , 


■ / 1/7 1 \ di 

Py = esm(/)i -h - «i) + 7^ 

'‘B 2 ‘■B 3 

Py = €sm4>5 - j^{d2 - di) - 

‘‘Bi ‘■B 3 


(43) 

(44) 


It is clear that the refraction angles (j)^ at the interfaces are obtained by requiring conservation of the 
momentum py. This leads to a simplified expression of these angles in terms of 4>i 


1 ■ 1 ^ / 7 7 \ 

sm05 = sm0i -b 3^(02 - di) -b 




^^B3 


and therefore we characterize our waves by introducing a critical angle (j)c 


= sm 


-1 


l-b2cii 



(45) 


(46) 


This tells us that when the incident angle is less than (/>c, the modes become oscillating guided modes, 
while in the case when the incident angle is more than (pc, we obtain decaying or evanescent wave 
modes. 

In the forthcoming analysis, we will be interested in studying the situation where pi < pc- To 
simplify our task and proceed further, let us choose the interfaces separating regions as 


O-'Ua — 




^n-i I'Bi 


bri- — 


Sj\^ 


ijenj Ibj (e^ Ibj + SjPjh ,) 


(47) 


We match the wave functions at the boundaries {—d 2 , —di, di, d 2 ) as required by the first order nature 
of the Dirac equation. For this, we introduce the shorthand notations 

_ _ _.—d2 




= D 


(±V2( 


-d2 




+ kyls^ 


(48) 

(49) 
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the related symbols by letting —d 2 


-di, 


’'ta = 0((.,.,(±^ (-^ + 

)'-(»-.,)’)/'2 (="'^ ( 7 ^ + ‘^■''«»)) 

the related symbols ^^3 follow by letting —di —)■ di, 


(50) 

(51) 


(52) 

(53) 


the related symbols follow by letting di —)■ d 2 - Now, requiring the continuity of the spinor 

wavefunctions at each junction interface give rise to a set of equations which can be expressed in terms 
of 2 X 2 transfer matrices between different regions 




ttj+i 

bj+i 


(54) 


where Mjj+i is a transfer matrix that couple the wave function in the j-th region to the wave function 
in the {j + l)-th region. Finally, we obtain the full transfer matrix over the whole double barrier 
region, which can be expressed in an obvious notation as 


ai 

bi 


4 


“ n ^11+1 

i=i 





(55) 


The total transfer matrix M = M 12 ■ M 23 ■ M 34 • M 45 is a transfer matrix that couple the wave function 
in the incident region to the wave function in the transmission region. It can be expressed explicitly 
as 


M 


Mi2 


Tf23 


M34 


M45 


mil rnu \ 

m2i m22 J 

g ^Pxl ^2 ^'^Pxl ^2 


-1 


zie 


-ipxld2 _^*^ipx\d2 


®«2^1n2 ^n 2 Vln 2 

'^bn2Ctn2 “*^«2?ln2 


^n2d2n2 

®n2%n2 


( Ona^/lns 


'^bn2^2n2 

— ibn2^2n2 / 

] ' 

^ '^bnsCtns 

—ibnziinz 


^md2n3 


( <^n2Vtn4 

®n2^1n4 

'^bnsitns. 

— ibn3^2nz / 

] ' 

\ '^bn2Ctn4 

~ibn2iln4 


-1 


®n2%n4 ®fi2%n4 

'^bn2^2ni ~'^bn2^2n4 


oiPx5d2 


,-ipxS,d 2 


z^e 


ipxs,d2 _^*g-ipx5d2 


(56) 

(57) 

(58) 

(59) 

(60) 


We consider an electron propagating from left to right with energy elB 2 > then r = bi and t = a^, r and 
t being the reflection and transmission amplitudes, respectively. We have assumed an incident wave 



from left normalized to unit amplitude oi = 1 and 65 = 0 is the null amplitude due absence of left 
moving waves in transmission region. This will give rise to the following relations 


1 


t =-, r = -. (61) 

mil "111 

This last formulation will be much more adequate in dealing with periodic systems and applying Bloch 
theorem to find the associated energy bands. The above expressions can be written as 


t = 


1 

mil 




r = 


m2i 

mil 




(62) 


where (ft and ipr refers to the phase of the transmission and reflection amplitudes, respectively. After 
a lengthy but straightforward algebra, we can show that t in (62) takes the form 




A+X+ + A-x“ + i(A-x+ - A+X“) 


(x+)2 + ix-? 


where we have set 


(63) 


A’^ = {l + {qt)^-{Qif)sin{d2{pxi+Px5))-‘^qtqicos{d2iPxi+Px5)) 

A“ = (1 + (qf f - (gf )^) cos(d 2 (Pxi + Pxs)) + q{ sm.{d 2 {px\ + Pxb)) 

X"^ = -al^bl^/dnsD - al^bl^an^C + S2S3an2an3bn2bn3{qiBi + q^B 2 ) 

-ssal^ansbl^bmAiiq+q^ -q^qi) 

X ~ ~®n 2 ^n 3 /^fi 3 A' ~ ®n3^n2*^^3'^ “k S2S30n2Q'n3^n2^n3 (^1 Bi + q^ B 2 ) 

+S3al^an3bl^bn3{{qtqt - QIQ5)^^+ 

-^1 — ^712^^3^714 H“ /37l2 TtIs ^714 

-^2 — '7712^713^714 “1“ *^712^713/3714 

Bl = ^742 ('^^277137714 3“ ‘ 3742 * 3713 / 3714 ) 

B2 — 0-742 (*^712 *3713 *3744 “(“77127713*^714) 

c = fe+ 0712 *3714 (52^1 ^712/3712 “1“ O74277I2) '^2^712^1 /3712(^5 0712*^714 ^2^7127714) 

D — (l^ 0-742 0^744 (O742 0^742 ^2^712*3712) “(“ ^2^1 ^712*3712(^5 0742*^714 ^2^7127714) 

E — (07420^742 '® 2 ^ 7 i 2 * 37 i 2)(^5 0742^744 '52^7127714) ^ 2^5 0742*^744^]^ ^712*3712 

E = {^2Qi ^ 712/3712 3“ 074277i2)(^5 0742*^714 ^2&7i277i4) ^2^5 ^n2^n^bn2Q.\ Pn2 

Zj = (?+ + iq- 

Oiuj = Pirij dtrij ~ Vlrij ^ 2 nj 

^lnj^2nj ^lnj^2nj 

luj = Vluj^trij + 

^rij = P2nj^lnj + vtrij^lnj 
^rij = V2nj^2nj + 'd2nj^2nj- 


The phase shift can be expressed explicitly as 


(ft = arctan 


A-X+ - A+x~ 
A+X+ + A-x" 


(64) 
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with the quantities 


Ax'^-A+x = 2cos(/>i(x'^cos((/)i + (pa;i+pa;5)) - X sm{4>i - {p^i + p^^)) (65) 

A+x"^ + A“X“ = 2cos(/>i(x’^sin((?ii - {p^i + Pxb)) + X~ cos{4>i + {pxi +Px5))- (66) 


(67) 


Finally the transmission phase is given by 

^ -1 [x^ cos((/>i + {pxi + Pxb)) - X~ sin((/>i - {pxi + Pxb)) 

— tcin ——- — - - -r--- — - - -^ 

[X+ sin(())i - [pxi +Pxb))+X~ cos(0i + [pxi +Pxb))_ 

Now we are ready for the computation of the transmission T and reflection R coefficients. For this 
purpose, we introduce the associated current density J, which defines T and R as 


T = 


Jt 


tra 


R = 


Jref 


( 68 ) 


'^inc '^inc 

where Jinc, Jref and J^a stand for the incident, reflected and transmitted components of the current 
density, respectively. It is easy to show that the current density J reads as 

J = evpfjj^axfl^ (69) 

which gives the following results for the incident, reflected and transmitted components 


Jinc = 

Jref = euF(V'r)^^^V'r 

Jtra = ■ 


(70) 

(71) 

(72) 


The energy conservation 


Pxi + (Py + ~P 2 ~id-i — J 2 ) — 






Px5 + (Py + 72 “ (^2 — Ji) + j^diY 
‘■B 2 ''B 3 


(73) 


allows us to express the transmission and reflection probabilities in the following simple forms 

2 


Px5 


1 


T = 

Pxi |miiF' 

More explicitly the transmission coefficient T reads as 

4pa,5(cos0i)2 


R = 


m 2 i 


mil 


T = 


Pxi [(x +)2 + ix-f] 


r,4 lA „2 .2 J.2 x2 x2 

— \ 2 ] ^ 712 ^ 712 ^^ 3 ^ 723 ^ 712 ^ 713^724 


(74) 


(75) 


Obviously, R and T are not independent, they are related through the unitarity requirement T+i2 = 1 
that is clearly shown in Figure 2a. Note that (45) implies that for certain incidence angles cjii the 
transmission is not allowed. In fact for 

Jbz < 7 — id 2 -di) + (76) 

1^2 ^B 2 \ ^-63 / 

all waves are completely reflected. 

We show the numerical results for the transmission, reflection coefficients and the GHL shifts in 
Figures 2, 3, 4, 5, 6 , for several parameter values (e, V 2 , U 3 , pj , di, ^ 2 )- For instance a typical value 
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of the magnetic field, say B 2 = dT, the magnetic length 13/^2 = 13nm, and = 1 corresponding to 
the energy E = 44mey [7] , these typical values will serve to normalize the various variables. The polar 
graph, Figure 2b, shows the transmission as a function of the incidence angle, the outermost circle 
corresponds to full transmission, T = 1, while the origin of this plot represents zero transmission. 
Requiring that eZsj = 3.7, d 2 = di, Ib 2 = ^2 = V 3 = 0, ^ = {0.5,1.5,3,3.67} and = 0 

reproduces exactly the result obtained in previous work [7]. Similarly, the transmission as a function 
of energy e for fixed -j^ = 0.8, ^ = 0.2 and = 0.6, i.e. ^ (^) ~ 1-156, shows that 

the transmission vanishes for — 1-156. 



(b) 


V3/s, 



Figure 2 : (a): Graphs depicting the reflection R (green line) and transmission T (red line) coefficients 
as function of energy potential v^Ib^ foi' the monolayer graphene barriers with = 0.2, ^ = 0.8, 

V2IB2 = 26 , elB2 = 30 , kylB2 = 1 ' = 2 and ixIb2 = 4 . (b): Polar plot of a curve with radius 


(transmission T) as a function of angle with = {15,7,5,1.35}, = 0.2, = 0.8, j 

V2IB2 ~ 3.1, V3IB2 — 1-2 and ^Ib 2 ~ 2. 

4 GHL shifts for double barriers 


J!3 _ 


= 0 . 6 , 


In this section, we shall turn to the Goos-Hanchen like (GHL) shifts in graphene by considering an 
incident, reflected and transmitted beams around a given transverse wave vector ky = ky^ and angle 
of incidence (j)i{kyQ) G [0, |], denoted by the subscript 0. These can be expressed in integral form 


1 


/ + 00 

dky f{ky - ky,) ^ + 

-00 

/ + 00 / 1 

^ dky r{ky) f{ky - ky,) i 


(77) 

(78) 


The reflection amplitude can be written as r{ky) = because of the ^-component of wavevector 

kxi as well as (pi are function of ky, where each spinor plane wave is a solution of (4). The angular 
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spectral distribution f{ky — ky^) can be assumed of Gaussian shape 


fiky 


ky,) = Wye-^y^^y-’‘yo^" 


( 79 ) 


where Wy being the half beam width at waist [18]. We can approximate the fey-dependent terms by a 
Taylor expansion around feyg and retain only the first order term to obtain 


Mh) - Mkyo) + . (ky- ky^) 


kxi{ky) ^ kxi{kyQ^ 


dky 

dk 


-HO 


xl 


dk. 


(ky kyg). 


-HO 


(80) 

(81) 


The transmitted wave takes the form 


^t{x,y) 


f' + OO 


dky t{ky^ f {.^y 


kyo) 


1 

J4>5{ky) 


(82) 


where the transmission amplitude t{ky) = is calculated through the use of boundary conditions. 

In order to determine the GHL shifts of the transmitted beam through the graphene double barriers, 
we adopt the following definition [21,22] 


Q - 

* dky 


KyO 


(83) 


s, 




Figure 3: The GHL shifts and the transmission as function of energy eZsj for the monolayer graphene 
barriers with ^ = 0.8, V 2 IB 2 = 30, V 3 IB 2 = 30, fey/sj = 1, = 0.5, ^ = 0.7 where {y 2 lB 2 = O' 

k 3 lB 2 = 8) color red, (/i 2 ^B 2 = 0, = 0) color green and (^ 2^52 = 8, = 0) color blue. 
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In Figure 3, the above transmission and GHL shifts are shown versus energy for different 
parameters of our system = 0.8, ^ 2^53 = ^ 0 ,V 31 b 2 = ^ 0 ,kylB 2 = Ij = 2, ^ = 0.7^ with zero- 
gap (/U 2 /S 2 = = 0): green color and finite gap {H 2 IB 2 = 0) color and {h 2 ^b 2 = 8, 

/^ 3 ^B 2 =d): blue color. It is clearly seen that GHL shifts are oscillating between negative and positive 
values around the critical point eIb = V 2 IB = The quantity kylB = m* plays a very important 

role in the transmission of Dirac fermions via obstacles created by a series of scattering potentials, 
because it is associated with the effective mass of the particle and hence determines the threshold for 
allowed energies. However, in the presence of an inhomogeneous magnetic field in the regions |x| < d 2 , 

it reduces this effective mass to i Py + 7^ (di — ^2) — ) in the incidence region while it increases 

V B2 'S3 J 


it to { Py — Tr~ {di — ^ 2 ) + 72 ^ ) in the transmission region. The allowed energies are then determined 

V B2 'S3 J 

by the greater effective mass condition: 


e^i32 > Py 


1 , , , , di 

Y2r\di — ^ 2 ) + 

‘‘B2 ‘Bs 


s, 




Figure 4 : The GHL shifts and the transmission as function of the potential V3IB2 for the monolayer 
graphene barriers with ^ = 0 . 8 , V2IB2 = 27 , = 30 , kylB2 ~ ~ ~ 0.78 where 

{P2^B2 = 0, P3IB2 = 8) color red, {p2^b2 = 0, P3IB2 = 0) color green and {P2IB2 = 8, P3IB2 = 0) color 
blue. 

The above GHL shifts and transmission are plotted in Figure 4 in terms of the potential V 3 IB 2 for 
some values of the physical parameters. It is clearly seen that St is oscillating between negative and 
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positive values around the critical point V 3 IB 2 = such points the transmission vanishes for 

€Ib 2 — — '*^3^52 < e^-B 2 + f^ 3 ^B 2 oscillates otherwise. 



Figure 5: The GHL shifts and the transmission as function of energy for the monolayer graphene 
barriers. (a)/(b) with {V2IB2 = 30, V3IB2 = 15 )/( u 2 /_ b 2 = V3IB2 = 30), with ^ = 0.8, kylB2 = 1. 

^ = 2, fJ, 2 lB 2 = fJ' 3 ^B 2 ~ = 0.19 (blue line), ^ = 0.4 (green line) and ^ = 0.7 (red line). 

In Figure 5, the transmission and GHL shifts are shown versus energy e/^j- can notice that 
at the Dirac points {elB 2 = ^ 2 ^B 2 5 ^^B 2 = 'V 3 IB 2 ), the GHL shifts change their sign. This change shows 
clearly that they are strongly dependent on the barrier heights. We also observe that the GHL shifts 
are positive as long as the energy satisfies the condition > V 2 IB 2 > V 3 IB 2 (Figure 5a) and negative 
for elB 2 < V 2 IB 3 < V 2 IB 2 (Figure 5b.) 

In Figure 6 , we analyze the transmission coefficients versus the potential V 3 IB 2 and V 21 b 2 - Ici do¬ 
ing so, we fix the energy elB 2 = 30 and choose a value of then we compute the transmission as 
shown in Figure 6 a. We notice that the transmission decreases if ^ increases and then vanishes while 
Figure 6 b shows different behavior. Note that, the Dirac points represent the zero modes for Dirac 
operator [16] and lead to the emergence of new Dirac points, which has been discussed in different 
works [23,24]. Such points separate the two regions of positive and negative refraction. In cases where 
V 2 IB 2 < c^B 2 V 2 IB 2 > c^B 2 (respectively V 3 IB 2 < e^S 2 ^-^Ld V 2 IB 2 > ^^B 2 )j the shifts are respectively 
in the forward and backward directions, due to the fact that the signs of the group velocity are opposite. 


14 













































































Figure 6: The GHL shifts and the transmission as function of energy potential V 3 IB 2 and V 2 IB 2 foi" the 
monolayer graphene barriers, (a): elB 2 = 30, V 2 IB 2 ~ ^ ~ ^y^B 2 = 1' = /^ 3 ^B 2 = 

4, = 0.2 (blue line), ^ = 0.5 (green line) and ^ = 0.78 (red line), (b): = 30, V 3 IB 2 = 32, 

^ = 0.8, kylB 2 = 1' = 2, plIb 2 ~ ~ (blue line), ^ = 0.12 (green line) and ^ = 0.02 

(red line). 


5 Conclusion 

To conclude, we have studied the transport of electrons in graphene scattered by double barrier in the 
presence of an inhomogeneous magnetic field. We obtained the solutions for the energy spectrum taking 
into account the conservation energy and noticed that for certain incidence angles the transmission is 
not allowed for e/s, < Py—^[di—d 2 ) + ^d\. However, the transmission probability T does not vanish 

‘S2 ‘S3 

in general, we also found that, in contrast to electrostatic barriers, magnetic barriers are able to confine 
Dirac fermions. This allowed us to calculate the GHL shifts of reflected and transmitted electron beams 
in a graphene double barrier structure in the presence of an inhomogeneous magnetic field. We also 
established some correlation between the electronic transport properties of Dirac fermions with the 
GHL shifts. 

The numerical data showed how these shifts behave in relation to the transmission probability T. 
It is found that the GHL shifts can be modulated by the incident energy potential energies V 2 IB 2 
and V 31 b 2 - The GHL shifts still change sign, but the point where it changes sign has been displaced 
to the left and the absolute value of the maximum of the shifts increased as well. Thus we seen that 
the GHL shifts in the transmission region can be either negative or positive. 
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